

for lnd = 1:length(beh)
    for knd = 1:length(beh(lnd).M)
        for ind = 1:length(beh(lnd).M(knd).data(10,:))
            if beh(lnd).M(knd).data(10,ind)>2.5
               beh(lnd).M(knd).angle(ind)=(beh(lnd).M(knd).data(10,ind)-2.5)/2.5*180;
            end
            if beh(lnd).M(knd).data(10,ind)<2.5
               beh(lnd).M(knd).angle(ind)=(abs(beh(lnd).M(knd).data(10,ind)-2.5))/2.5*(-180);
            end
        end
    end
end

figure,
id = 0;
for ld = 1:length(beh)
    for kd = 1:length(beh(ld).M)
        id = id+1;
        subplot(length(beh),length(beh(ld).M),id), hist((beh(ld).M(kd).angle),100)
        xlim([-185 185])
        ylim([0 400000])
    end
end

for lnd = 1:length(beh)
    for knd = 1:length(beh(lnd).M)
        beh(lnd).M(knd).polar = beh(lnd).M(knd).angle * pi/180;
    end
end

for lnd = 1:length(beh)
    figure, rose(beh(lnd).M(3).polar,1000)
end


h = [subplot(231) subplot(232) subplot(233) subplot(234) subplot(235) subplot(236)] 
axis([h], xlim)
